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RESEARCH MEMORAUDIM 


LONGITUDIML STABILITY AND CONTROL CHARACTERISTICS OF A 
CANARD MISSILE CONFIGURATION FOR MACH NTMBERS FR(M 
1.1 TO 1.93 AS DETERMINED FRCM FREE-FLICTT 
AND WIND-TUNNEL INVESTIGATIONS 
By Howard J. Curfman^ Jr.^ and Carl E. Grigsby 

SUMMARY 


The supersonic longitudinal stability and control characteristics 
of a canard missile configuration as determined from flight tests of 
rocket -powered models and supersonic wind-tunnel tests have been com- 
pared. The flight tests covered a Mach number range from 1.1 to about 
1.6 representing a Reynolds number range of about 10 x 10° to 20 x 10°, 
based on wing mean aerodynamic chord. The wind-tunnel tests were at 
Mach numbers 1.62 and 1.93 with Reynolds numbers of O .65 x 10° and 
0.58 X 10°, respectively. 

The results from the two techniques showed good correlation for 
the slopes of the lift and moment coefficients against angle of attack 
and for aerodynamic- center location. The nonlinearities in these 
stability derivatives were present to the highest Mach number attained . 
The control effectiveness of the deflected canard fins was shown to be 
practically a pure moment with very little change in lift, the over -all 
lift effectiveness of the configuration "being due to the ability of the 
canard fins to change the angle of attack. The average moment effec- 
tiveness as Indicated by the flight tests was higher than that shown by 
the wind-tunnel results . The use of the stability derivatives as 
determined by both techniques in predicting the longitudinal frequency- 
response characteristics gave satisfactory agreement . 


INTRODUCTION 


The aerodynamic characteristics of missile configurations have 
received intensified study in recent years. This paper compares 
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cheiracteristics of a canaxd missile configioxation as determined by free- 
fli^t tests using rocket -powered models and supersonic wind-tunnel 
results . Eie wind-tunnel data are presented herein, and the flight results 
are as presented in references 1 and 2. This configuration has been used 
as a research vehicle in the transonic and supersonic speed regions, and, 
in addition to obtaining aerodynamic information, evaluations of automatic 
stabilization systems have also been conducted. References 1, 2, and 3 
give information on the longitudinal and roll characteristics of the 
configuration and include the control effectiveness of. canard elevators 
and wing -tip ailerons. In the material presented herein, only those longi- 
tudlnaJ. data from both techniques which can be compared directly are given. 


STfMBOLS , 

lift coefficient, Lift/qS 

Cjij pitching-moment coefficient. Moment /qSc 

Cp drag coefficient, Drag/qS 

q dynamic pressure, Ib/sq ft 

5 total wing area in one plane, formed by extending leading 

and trailing edges to center line, sq ft- 

V velocity, ft/sec 

c wing mean aerodynamic chord, ft 

a angle of attack, deg 

6 canard-fin deflection 

M Mach number 
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m 


Sc_ 


•"nin 


exposed canard -fin area in one plane, sq ft « 
moment of inertia in pitch,' slug-ft^ 
mass, slugs 


m 


2V 


p m 


a 


OJ 


|9/8| 


2V 


pitching velocity, radians/sec 

time rate of change of angle of attack, radians/sec 
forcing frequency, radians/sec 
amplitude ratio, radians per sec/radian 


DESCRIPTION OP MODELS 


The flight model of reference 1, designated model A herein, had 
cruciform wings and all -movable canard control surfaces of approximately 
60° delta plan form. The aerodynamic arrangement was such that the 
horizontal wings and controls were in the same plane . These surfaces 
were mounted on a cylindrical body of fineness ratio 16.3 having ogival 
nose and tail sections of about 3*5 calibers. The wings and canard 
surfaces had modified double -wedge airfoil sections with constant thick- 
nesses corresponding to a thickness ratio of about 3 percent at the wing- 
body juncture and about 3 percent at the canard root chord. 

The flight model of reference 2, designated herein as model B, was 
like model A in many general aspects but was not geometrically similar. 
Model B was approximately 0.875 scale of model A. The wind-tunnel model, 
model C, was more nearly geometrically similar to model A, differing 
slightly in the length of the cylindrical portion of the fuselage. The 
tunnel model was approximately O.0875 scale of model A. To illustrate 
more clearly the differences between the three models, figure 1 and 
table 1 are prepared; in these presentations the ratio of diameters is 
the linear scale factor used to convert dimensions of models B and C to 
correspond to those of mdSel A. The boattarlling of models A and B was 
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impractical la the -wind tunnel, and in figure 1 the dotted lines indi- 
cate the relative sting-windshield position and size . Figure 2 shows 
the details of the exposed wing and canard surfaces with dimensions 
corresponding to model A. 


TEST TECHNIQUES 
Rocket -Powered Model Tests 


In the rocket-powered model tests of references 1 and 2 , the model 
contained no sustainer or internal rocket motor hut was boosted to the 
maximum supersonic velocities by means of solid -propellant rocket motors. 
During the coasting fli^t from supersonic to subsonic Mach numbers 
following booster separation, data were continuously transmitted by 
telemeter and recorded at ground stations- In general, these data 
included longitudinal and normal accelerations, angle of attack, and 
total and static pressures. In addition, radar tracking data and radio- 
sonde infonnation were available. 

The longitudinal stability and control derivatives were obtained by 
analyzing the model response to programmed step inputs of the canard longi- 
tudinal control surfaces. The period and dsimping characteristics of the 
transients and the measxarements of steady state or trim values allowed a 
determination of the stability and control derivatives from the use of the 
general equations of motion considering two degrees of freedom longi- 
tudinally for each available transient. A detailed discussion of this 
technique is presented in the appendix of reference 2. 

For model A the programmed control-surface deflections were of a 
square wave with deflections of 2.68° and -4.42°. The reason for this 
approach was that such a method would generally give trim values of angle 
of attack of sli^tly different magnitudes and could thus indicate the 
presence, if any, of nonlinear regions in the various parameters. For 
model B the square-wave input involved control deflections of 


Wind-Timnel Tests 

All the wind-tunnel test resiilts presented herein were obtained in 
the Langley 9-inch supersonic tunnel which is a continuous-operation 
closed -circuit type in which the stream pressure, temperature, and 
humidity conditions can be contro3J.ed. Different test Mach numbers are 
provided by Interchanging nozzle blocks which form test sections 
approximately 9 inches square. Throughout the wind-tunnel tests, the 
moisture content in the tunnel was kept sufficiently low so that the 
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effects of condensation in the supersonic nozzle were negligible. 

Eleven fine-mesh turbulence -damping screens are provided in the 
relatively large area settling chamber just ahead of the supersonic 
nozzle. 

The wind-tunnel model was sting-supported and the sting and sting- 
windshield arrangement used in these test is illustrated in figure 3* 

At each angle of attack^ the model, sting, and sting windshield were 
translated across the tunnel so that a fixed point on the model (center- 
of-gravity location of 0.586c ahead of leading edge of c) could be kept 
on the center line of the tunnel. By use of this arrangement, configu- 
rations, which at an angle of attack of 0° were free from shock 
reflections, could be tested through the angle range of il5°. Through- 
out the tests the gap between the rear of the movable windshield was 
less than 0.012 inch. The sting and sting -wind shield arrangement, 
illustrated in figure 3^ tms been shown by other tests to have negligible 
effect on the aerodynamic characteristics of the model. 

Measurements of lift, drag, and pitching moment were made by means 
of external self-balancing mechanical scales through an angle -of -attack 
range of -5° to 15° for canard-fin deflection angles of 0°, 3 °} 6°, 10°, 
and 15° . Measurements of the pressure in the sting-windshield and balance 
enclosing box, which tests have shown to be equal to the model base 
pressure, were made and the drag results were correcJ;ed to the condition 
of base pressure equal to, stream pressure. An optical system employing 
a small mirror mounted in the rear of the model was used to measxire 
angles of attack. * 


Test Conditions 

The Reynolds number of the flight tests of model A (ref. 1), based 

on the wing mean aerodynamic chord, varied from 10 x 10^ to 20 x 10^ for 
the Mach numbers from 1.1 to 1.63, respectively. For model B (ref. 2) 

the Reynolds number variation was from 10 x 10^ to 15 x 10^ for Mach 
numbers from 1.1 to 1.^5, respectively. The Reynolds numbers of the 

wind-tunnel tests at Mach numbers of 1.62 and 1.93 were 0.65 X 10^ and 
0.58 X 106, respectively. 


Accuracy 

The accuracy which may be expected from the free -flight tests of 
rocket -powered models is affected by many considerations, including 
telemetering and recording errors, velocity errors resulting from both 
radar and pressure measurements, reading errors, and others. 

# 
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Reference 2 ^contains a discussion on accuracy and gives the following 
as the estimated prohahle errors in the stability derivatives and 
coefficients of interest in the present paper for M = 1.40: 


Percent of given value for - 


Si 



^D 

2 

4 

ho 

6 

h 


^'Thile no discussion of the accuracy of the flight tests results of 
reference 1 is detailed therein, the accuracy is believed' to be reasonably 
the same as that given for the flight tests of reference 2, since the 
technique of data reduction and instrumentation was essentially the same. 

The precision of the wind-timnel resiilts has been evaluated by 
estimating the uncertainties in the balance measurements Involved in a 
given quantity and combining these errors by a method based on the theory 
of least squares. 

A summary table of precision estimates follows : 


Lift coefficient, Cj^ *0.0002 

Drag coefficient, *0.0003 

Pitching-moment coefficient. Cm , ±0.001 

Angle of attack, a, deg . . ±0.01 

Canard -fin deflection angle, 6, deg ±0.25 


RESULTS AND DISCUSSION 


The flight results presented herein were taken directly from 
references 1 and 2. The wind-tunnel results at Mach numbers 1.62 eind 1.93 
are presented in figure 4 with the pitching-moment coefficient referenced 
to a center -of -gravity location of 0.568c ahead of the leading edge of 
the mean aerodynamic chord. The various characteristics of the curves 
of figure 4 are discussed in detail in the sections which follow. 


Lift 


The variation of the lift-curve slope Cj^ with Mach number is 


a 


shown in figure 5- The wind-tunnel points are given for 5 = 3° ±he 
small angle-of -attack region (0< a< 2°) and for 6 = 6° in the angle* 
of -attack region greater than .2°^ to correspond more closely to flight 
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results . The correlation of flight and wind-tunnel results for 

is considered quite good. The wind-tunnel data confirm the nonlinearity 
in lift-curve slope indicated hy the flight test of reference 1 and 
reveal a slope change of about the same magnitude. The tvinnel results 
indicate that this nonlinearity persists at the highest Mach number 
tested. The wind-tunnel curves of figure 4 show that the nonlinear 
lift -coefficient curve could in general he represented hy two linear 
segments with the break at about a = 2°. Similar nonlinear character- 
istics have been observed on other canard configurations (see^ for example, 
ref. 4) . The reduction in lift-curve slope with Mach number was nearly 
Identical with that calculated for a thin triangular wing from linearized 
theory. 


Static Stability 

The variation of aerodyuamic -center location with Mach number as 
determined from the flight and tunnel tests is shown in figure 6. The 
wind-tunnel points are given for 6=0°, 3°^ and 6° and were taken at 
trim from the slopes of the Cjjj - curves. It is noted that where 

the trim a is near zero, the aerodynamic -center location is somewhat 
forward from the positions indicated for the trim values of a for the 
larger, control deflection; this is due to the downwash effects being 
greater at angles of attack very near zero. Since the downwash effects 
are dependent on the angle of attack, it is expected, of course, that 
the indicated aerodynamic -center location as determined from plots for 
the two center-of-gravlty locations might vary at trim conditions for 
the same control deflection. In general, however, these differences 
amount to a maximum of about 0 . 03 ^ for 6=3° and 6° at trim for the 
two flight-test locations of the centers of gravity. Though the canard 
deflections of model B were -5°, the forward .center -of -gravity location 
resulted in smaller trim angles of attack than for model A, and the « 

results indicate a more forward aerodyneimic -center location. The data 
from model B, however, indicate a slightly greater forward shift than 
can be noted in the tunnel results . With this exception, just noted, 
the flight and wind-tunnel tests show good agreement. The forward 
movement of the aerodynamic center with Mach number was also found to 
be slightly more rapid than that indicated theoretically for the body— 
triangular-wing combination (no canard surfaces) discussed in reference 5* 

The static pitching-moment derivative is considered here in 

“a 

addition to the discussion of aerodynamic -center location, since the 
frequency of the short-period mode of motion of the configuration is 
almost entirely dependent on this parameter. The variation of this 
derivative with Mach number is presented in figure 7* The wind-tunnel 
data presented here were detennined for the corresponding center -of -gravity 
locations with respect to the mean aerodynemic chord as existed in the 
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flight tests , In the comparisons shown, the wind-tunnel data were used 
for 6 = 0°, 3°^ 8Jid 6° with the slopes at. trim presented, and . the 
agreement is quite reasonable. 

The nonlinearity in indicated hy the flight tests of model A 

and hy the tunnel tests is also noted for the canard configurations of 
reference 4 and is clearly shown from the wind-tunnel results in figure h. 
This characteristic is more pronounced in the regions of small angle of 
attack. The nonlinearity of Cm against a has been noted on many 

canard configurations with inline surfaces and is the well-known result 
of the downwash field at the rear wings. The reduction in stability 
noted in the tunnel tests at angles of attack greater than 10° is 
probably due to sepaxation effects. The effects of similar changes in 
Cjj^ on the problem of applying automatic control equipment to this 

configuration have been studied in reference 6, and the considerations 
therein lead to the conclusion that the magnitude of the slope changes 
in the region of small angles of attack as shown herein should cause no 
great problem in the automatic stabilization and control of this 
aircraft . 


Control Effectiveness 

The flight tests of references 1 and 2 showed that the horizontal 
canard surfaces produced very small changes in total lift but were 
effective in producing a pitching moment to change the angle of attack. 
These results are confirmed by the wind-tunnel tests . The very low 
values of Ct are believed to be due to the effect of the downwash 

from the deflected canard surfaces on the wings, the change in lift on 
the wings due to doimwash being opposite and nearly Identical to the 
lift on the canard fins . Although Cj^^ was the least accurate of the 

parameters presented in reference 2, there was reported therein 
negative values of this derivative in the supersonic range to M = 1.44. 
The effect of the downwash from the deflected canard fins onto the wings 
also results in an increase in the pitching moment since the wings and 
canards are on opposite sides of the center of gravity. Hence, the 
effectiveness of the deflected canard fins can be described as producing 
practically a pure moment or couple while of Itself causing very little 
change in lift, the over -all lift effectiveness of the configuration 
being due to the ability of the canards to change- the angle of attack. 

The wind-tunnel results of against 6 at constant a are 

presented in figure 8 for both test Mach numbers. The slopes of these 
curves are noted to be quite small; this plot indicates, however, that 
the slope trend is positive, the largest value being about 0.001. Values 
of Cl from -0.0026 at M = 1.1 to -0.0013 at M = 1.44 are shown in 



reference 2. Because of the accuracy Involved in determining this 
parameter by both techniques, perhaps the lift effectiveness of the 
canards can be summarized as being negligible in magnitude but with 
indications that the slope Ct_ may be positive or negative for the 


configuration with the present canard-wing locations and relative areas 
and spans of canards to wings. References 8, and 9 are results of 
supersonic wind-tunnel tests' of other canard configurations ; these 
references show small values of each case less than 0.003. It 

might well be mentioned that the wind-tunnel tests referred to in the 
preceding discussion were all conducted at a Reynolds number of less 
than 1 X lo6. Since the flight tests of reference 2 were in the 
Reynolds number range of from 10 x 10^ to 15 X 106, and if the accuracy 
involved does not affect the slope trend, the possibility exists that 
the negative values of reported therein (and also indicated in other 

unpublished flight-test data) are a result of scale effects, but are 
still very small in magnitude in the low supersonic range considered 
here . 


The wind-timnel results of Cm against 6 at constant a are 
presented in figure 9 and are cross plots from faired curves. These 
curves indicate changes in the slope Cm^ with canard deflection, 

.angle of attack, and Mach number. In general, there seems ' to be a 
slight reduction in Cmg at the large canard deflections and at the 

higher angles of attack for all canard deflections. The curves also 
are fairly smooth for M = 1.62 but tend to be more nonlinear at the 
higher angles of attack at M = 1.93 ♦ 


The curves of figure 9 also show the reduction, in the control 
effectiveness parameter as the Mach number is Increased. Figure 10 

presents Cmg against Mach number. The data for models A and- B are 


also shown as presented in references 1 and 2, respectively. The wind- 
tunnel results are shown for both center -of -gravity locations at angles 
of attack of 2° and 4° and are faired values of C_ for the lower 

mg 

deflection range. While the tunnel results tend, in general, to appear 
somewhat low compared to the flight results, the over -all picture of the 
reduction in Cjjj with Mach number is apparent. It is emphasized that 


the flight results were obtained from changes in trim conditions and the 

measured values of the static stability derivative Cm ; hence, for the 

cc» 

noted nonlinear variations of Cm with a, this method of linear 
analysis can only afford a mean or average value of Cm„ • Thus, if it 

mg 


is only assumed that the flight data are indicative of the mean or 
average control effectiveness, figure 10 shows that the reduction in 
control effectiveness with Mach number noted from the flight tests up to 
M = 1.6 is continued to the.. higher. Manh-ujunib ftrj 1.93* 
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In figure 11 tlie wind-tunnel results of against 5 are 

presented for the two center-of -gravity locations. These curves indicate 

that the rate of change of trim a with canard deflection becomes 

appreciably reduced as the control deflection is increased. No flight 

results are presented for comparison because of the flight methods of 

obtaining comparable data. In the flight tests a, . against Mach 

trim 

nimiber was found for the two values of control deflection and the ratio 
^trlm 

= was formed as a function of Mach number. The known presence of 

nonlinearities in against a and model out-of-trim precludes the, 

use of this ratio as being directly comparable to the slope measured in 
figure 11. It is apparent, therefore, that the flight technique cannot 
show the true shape of the against 6 curve. 


Application of Results 

Since the combination of the aircraft -response characteristics 
with those of automatic -control components is a necessity in msmy current 
applications, a comparison of the frequency-response curves of the air- 
craft as predicted by derivatives obtained from fli^t and wind-tunnel 
tests is of interest. Consideration is given herein to the ratio of 
pitching-velocity response to canard elevator deflection ©/& as 
determined from equations involving two degrees of freedom longitudinally. 
This expression in terms of the differential operator D = d/dt is 
(see ref. 10) • 

i = K(tD + 1) 

^ + 2t(D J) + CD ^ 

^ n n 

where 


- Vi'S 

^y mV 
qSc qS 


P mV 


- Ct Cm — 
■^6 a 2V 


T 
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2H = 



c 

2V 


2 ^ -Clu 2V ^ 

•^y mV -*T 

qSC qS qS^ 


where o)^ and ^ correspond to the conventional undamped natural 

frequency and the damping ratio, respectively, and all stability 
derivatives are in radian measure. 


Assiaming Cj;^ to be zero, it is noted in the above expressions 
that T becomes independent of Cjjj leaving only the gain factor K 

, 5 

being changed by the control effectiveness j therefore, the differences 


in 


Cjjjg as indicated from flight and tunnel tests produce only a scale 


change in the frequency-response function and do not affect the shapes 
of the frequency-response curves. 


The calcxilations were made for a Mach number of 1.6 at sea level 
using the mass and inertia values given in reference 1, namely, 
m = 5-05 slugs and ly = 31 -3 slug-feet2. Stability derivatives were 
defined for small and large angle-of-attack regions and may be read from 
figxires 5 and 7* The damping-in- pitch derivatives were those presented in 
reference 1 for the two angle-of-attack regions. These frequency-response 
curves are presented as amplitude and phase against cn, the forcing 
frequency, in figure 12. It is noted that the K factor has not been 
included in these resxilts. 


These curves can illustrate effects of the nonlinearities noted in 
the stability derivatives, as well as the magnitude differences shown 
from flight and wind-tunnel data. The effects of the nonlinearities are 
compared by the curves for high and low, angles of attack and illizstrate 
that the location and magnitude of the resonant peaks are slightly 
different. The difference in resonant frequencies noted for the high- 
angle-of-attack comparisons is accounted for by the differences in 
values of Cj^; the major portion of the peak -magnitude variations is 


caused by the damping factors. The maximum phase difference occurs near 
resonance. For many purposes of analysis these results are adequate. 


;v 
1 . 
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CONCLUDING REMAEKS 


A comparison has been made between the supersonic longitudinal 
stability and control characteristics of a canard missile configuration 
as determined from free-flight tests of rocket -powered models and tests 
in a supersonic wind tunnel. The flight-test range included Mach niun- 
bers from 1.1 to about 1.6 representing a Eeynolds number range from 

about 10 X 10^ to 20 x lO^j based on wing mean aerodynamic chord. The 
wind-tunnel tests were at Mach numbers of 1.62 and 1.93 with Eeynolds 

numbers of 0.65 X 10^ and 0.58 x 10^, respectively. 

Good correlation was shown between the two techniques for slopes 
of lift and moment coefficients against angle of attack and for 
aerodynamic -center location. The nonlinearities in these static 
stability derivatives due to angle of attack noted in the flight tests 
were also shown in the wind-tunnel tests and were present at the highest 
Mach number attained. 

The control effectiveness of the deflected canard fins can be 
described as producing practically a pure moment or coupld while of 
itself causing very little change in lift, the over-all lift effective- 
ness of the configuration being due to the ability of the canard fins 
to change the angle of attack. The average moment effectiveness of the 
canard fins as indicated by flight tests tended to be higher than shown . 
by the wind-tunnel results. 

The use of stability derivatives determined from flight and tunnel 
tests in predicting longitudinal frequency-response characteristics showed 
satisfactory agreement and would be adequate for many purposes of analysis. 

Langley Aeronautical Laboratory 

National Advisory Committee for Aeronautics 
Langley Field, Va. 
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TABLE I.- COMPAEISON OF FLIGHT AMD WIND-TUNMEL MODELS 


Characteristics 

Model A 
( ref . l) 

Model B^ 
(ref. 2) 

Model C®' 
(wind tunnel) 

Diameter, in 

8 

8 

8 

c , f t 

1.766 

1.703 

1.766 

S, sq ft 

4.110 

3.780 

4.110 

Sf, sq ft 

0.282 

0.251 

0.282 

Span, in , . 

37-000 

35'4 oo 

37.000 

Length, in 

130.375 

130.510 


Span/Diameter 

4.625 

4.430 

4.625 

c.g. location, 

c ahead of L.E. c 

0.384 

0.622 

0.568 

Length from trailing edge 
of wing to trailing edge 
of canard fin, in. . . . 

63.875 

64.230 

62.875 


^Models B and C dimensions have been scaled to model A by 
the ratio of diameters. 
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Figure 1.- Compariaon of geometry of fli^t models (refs. 1 and 2) and 
wind-tunnel model. Models B and C are drawn scaled to the diameter 
of model A. All dimensions are In Inches. 
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(a) M=1.62. 

Figure 4.- ¥1116.-111111161 resulls. Pllclilng 111011161118 are aloul a cenler of 
giravlly localed 0. 568?T ahead of Ihe leading edge of c. 
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Figure 12.- Comparison of flight and wind-tunnel results in predicting 
the pitching velocity per imit canard-fin-deflection frequency 
response. 
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